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INTRODUCTION 


This  introductory  discussion  has  much  in  common  with  our  previous  annual  report  of 
November  6,  1996,  since  the  work  has  continued  with  the  same  goals.  We  have  requested 
a  no-cost  extension  for  12  additional  months  to  complete  the  work. 

In  its  present  form,  magnetic  resonance  imaging  (MRI)  produces  images  by  mapping 
the  hydrogen  nuclei  in  the  tissues  of  the  body.  However,  certain  portions  of  the  body, 
most  notably  the  lungs,  have  remained  difficult  to  image  using  conventional  MRI.  Infor¬ 
mation  about  lung  function  has  been  even  more  difficult  to  obtain  with  presently  available 
radiological  techniques.  A  new  implementation  of  MRI  using  laser-polarized  noble  gases 
has  recently  been  demonstrated^’^  wherein  lasers  are  used  to  enhance  the  MR  signal  from 
noble  gases  such  as  ®He  and  ^^®Xe,  making  them  easily  observable  in  a  conventional  MRI 
scanner.  Initial  experiments  have  already  yielded  spectacular  magnetic  resonance  images  of 
the  lungs  of  laboratory  animals  and  of  human  volunteers.  This  technology  should  provide 
functional  information  that  can  be  important  in  evaluating  and  treating  chronic  obstruc¬ 
tive  pulmonary  diseases,  pulmonary  embolisms,  emphysema,  asthma,  lung  cancer  and  a 
wide  variety  of  other  respiratory  problems. 

BODY 

Magnetic  resonance  imaging  with  laser-polarized  noble  gases  has  been  made  possible 
by  years  of  basic  physics  research  in  the  areas  of  optical  pumping  and  spin  exchange®, 
largely  supported  by  AFOSR  with  more  recent  assistance  by  DARPA.  Optical  pumping 
uses  circularly  polarized  light,  most  often  from  a  laser,  to  create  a  large  electron  spin 
polarization  in  a  vapor  of  rubidium  or  similar  atoms.  In  a  collision  with  a  noble  gas  atom, 
the  electron  spin  of  the  rubidium  atom  can  be  transferred  to  the  nuclear  spin  of  a  noble-gas 
atom.  Extremely  large  nuclear  polarizations  can  be  obtained  on  time  scales  ranging  from 
minutes  to  a  few  hours.  With  Titanium: Sapphire  lasers  and,  more  recently,  high-power 
AlGaAs  laser  diode  arrays,  it  is  possible  to  polarize  nuclei  in  substantial  quantities  of  gas 
(a  few  liters  at  STP).  This  is  enough  for  use  in  magnetic  resonance  imaging. 

In  magnetic  resonance  imaging,  an  image  is  reconstructed  from  the  radio  waves  pro¬ 
duced  by  processing  nuclear  magnetic  moments  fi.  For  images  of  human  beings,  the  pro¬ 
cessing  moments  are  almost  always  protons,  the  nuclei  of  hydrogen  atoms.  However,  the 
radio  signal  from  a  spin-up  nucleus  completely  cancels  the  radio  signal  from  a  spin-down 
nucleus.  The  observable  signal  in  magnetic  resonance  imaging  is  due  to  the  spin  polariza¬ 
tion  P  —  {Ni  —  Ni)/N,  that  is,  the  difference  between  the  number  of  spin-up  nuclei 
and  the  number  iVj  of  spin-down  nuclei,  divided  by  the  sum  N  =  +  Ni.  In  conventional 

MRI,  the  polarization  is  the  normal  Boltzmann  polarization,  P  =  fiB/kT,  equal  to  ratio 
of  the  magnetic  energy  difference  fiB  of  a  spin-up  and  a  spin- down  nucleus  in  a  magnetic 
field  B,  to  the  thermal  energy  kT.  In  practically  attainable  magnetic  fields  and  at  human 
body  temperature,  the  resulting  polarization  is  only  a  few  parts  per  million.  It  is  only 
because  of  the  huge  number  of  protons  in  most  organs  of  the  human  body  (mostly  proton- 
rich  water  and  fat)  that  the  radio  waves  from  this  part-per-miUion  proton  polarization  are 
strong  enough  to  create  an  image. 

By  contrast,  optical  pumping  and  spin  exchange  produces  noble  gas  nuclear  polariza¬ 
tions  of  order  1,  so  that  the  magnetic  resonance  signal  per  noble  gas  nucleus  is  five  orders 
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of  magnitude  larger  than  the  signal  per  proton  in  conventional  MRI.  Thus,  even  though 
the  density  of  nuclei  in  a  gas  is  one  thousand  times  smaller  than  the  density  of  protons  in 
tissue,  the  huge  increase  in  polarization  more  than  compensates  for  the  fewer  nuclei  and 
very  bright  images  of  laser-polarized  gas  can  be  obtained. 

The  first  gas  images,  of  the  excised  lungs  of  a  mouse,  were  made  in  1994  with  a  few  cc’s 
of  laser-polarized  ^^®Xe  by  a  team  of  researchers  from  Princeton  University  and  the  State 
University  of  New  York  at  Stony  Brook^.  Further  collaboration  between  Princeton  and 
researchers  at  Duke  University®  yielded  the  first  images  made  with  laser-polarized  ®He;  the 
first  lung  images  of  a  living  gmnea  pig  were  reported  early  in  1995.  A  team  from  Princeton 
and  Duke  Universities  produced  the  first  human  lung  image  at  the  Duke  Medical  School 
on  September  18, 1995.  Clinical  research  at  the  University  of  Virginia  Medical  School  and 
at  the  Duke  University  Medical  School  is  now  producing  several  human  lung  images  per 
day. 

MR  imaging  with  laser-polarized  noble  gases  ofiers  a  wide  variety  of  possible  clinical 
applications.  The  large  non-equilibrium  polarizations  allow  rapid  imaging  techniques  to 
be  used.  Presently,  a  single  “slice”  of  a  human  lung  can  be  imaged  in  much  less  than 
one  second,  so  real-time  imaging  of  lung  function  is  possible  for  the  first  time.  The  high 
resolution  offered  by  the  technique  and  the  intrinsic  three-dimensional  nature  of  MRI 
means  that  areas  of  compromised  lung  function  can  be  localized  to  a  much  higher  degree 
than  presently  possible.  Current  clinical  techniques  involve  the  use  of  inhaled  radioactive 
^®®Xe  (a  gamma  emitter)  to  produce  a  two-dimensional  projection  image  of  the  lung  with 
resolutions  on  the  order  of  1  cm  X  1cm.  The  technique  is  frequently  combined  with 
radioactive  ®®Tc,  injected  intravenously  and  imaged  in  a  similar  manner.  By  contrast, 
fully  three  dimensional  human  images  using  laser-polarized  ®He  are  routinely  obtained 
with  a  resolution  of  1  mm  X  1  mm,  and  much  better  resolution  has  been  demonstrated  in 
laboratory  animals. 

While  initial  experiments  have  concentrated  on  imaging  in  the  gas  phase,  laser- 
polarized  ^®®Xe  offers  the  further  possibility  of  imaging  the  blood  and  tissues  as  well. 
Xenon  is  highly  soluble  in  tissue  (especially  fatty  components)  and  will  dissolve  in  the 
blood  at  levels  of  order  20%  of  its  concentration  in  the  gas  phase.  The  polarization  of 
^®^Xe  survives  for  lO’s  of  seconds  in  the  blood,  so  imaging  of  ^®^Xe  in  almost  any  desired 
portion  of  the  body  will  likely  be  feasible.  Dissolved  ^®®Xe  has  a  precession  frequency  that 
is  easily  distinguished  from  the  gas-phase  frequency,  so  gas-phase  and  tissue  images  can  be 
recorded  separately.  Detection  of  pulmonary  emboli  (blood  clots  in  the  lung),  a  condition 
that  contributes  roughly  200,000  deaths  per  year  in  the  United  States,  is  envisioned  as  a 
procedure  which  could  make  excellent  use  of  ^®®Xe  imaging,  with  gas-phase  images  giving 
much  better  and  quicker  maps  of  lung  ventilation  than  conventional  radioactive  ^®®Xe, 
and  with  dissolved  ^®^Xe  signals  giving  better  and  quicker  information  than  radioactive 
^^Tc  on  blood  perfusion  of  the  lung.  Other  diagnostic  possibilities  are  certain  to  arise 
once  large-scale  clinical  research  commences.  Lung  function  is  an  important  concern  for 
the  DoD  in  many  areas,  including  the  lung  function  of  aviators  at  high  altitudes,  the  lung 
function  of  deep  sea  divers,  and  the  functioning  of  lungs  in  chemical- warfare  environments. 

Thus  far,  noble  gas  MRI  experiments  have  made  use  of  commercial  MRI  scanners 
which  have  been  retuned  to  the  noble  gas  precession  frequency  (a  relatively  straight-forward 
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and  inexpensive  task).  However,  because  the  signal-to- noise  ratio  of  the  MR  signed  from 
laser-polarized  noble  gases  is  independent  of  magnetic  field  strength,  imaging  with  much 
smaller  fields  will  be  possible.  A  dedicated  MRI  unit  for  noble  gas  imaging  would  be  small, 
inexpensive,  and  much  more  portable  than  most  conventioned  units. 

The  means  of  production  of  laser-polarized  noble  gases  are  far  from  optimized,  and 
research  is  continuing  to  learn  how  to  efficiently  and  cost-effectively  produce  the  large 
quantities  of  gas  needed  for  medical  and  other  imaging  applications.  Interactions  of  the 
noble  gas  atoms  with  surfaces,  as  well  as  further  subtleties  of  the  optical  pumping  and  spin 
exchange  processes  will  continue  to  be  explored  with  the  potential  for  substantial  gains  in 
production  capability. 

CONCLUSIONS 

The  aim  of  this,  work  as  outlined  in  the  Statement  of  Work  attached  as  Appendix 
A,  was  to  develop  as  rapidly  as  possible  the  capability  to  produce  large  amounts  of  nuclear 
spin-polarized  ^^®Xe  and  ®He  gas  for  magnetic  resonance  imaging  of  humans  and  for  other 
applications.  This  promising  new  diagnostic  method  uses  conventional  magnetic  resonance 
imaging  machines,  along  with  laser-polarized  ^^®Xe  and  ®He  gas,  to  make  images  of  human 
lungs  and  other  body  cavities.  Xenon  is  so  soluble  in  blood  and  tissue  that  it  will  be  possible 
to  image  other  organs  like  the  brain  with  laser-polarized  ^^®Xe.  The  research  has  been 
very  successful,  and  some  of  the  highlights  are: 

•  Many  ®He  images  of  human  lungs  have  now  been  obtained  at  the  Medical  Schools  of 
the  University  of  Virginia  and  Duke  University. 

•  Progress  has  been  made  in  understanding  the  operation  of  diode  laser  arrays  for  spin- 
exchange  optical  pumping  of  ^^^Xe  and  ®He  gas.  This  work,  has  been  facilitated  by 
the  development  of  a  powerful  new  imaging  diagnostic,  perfected  over  the  past  year, 
for  mapping  out  the  spatial  distribution  of  spin  polarization  of  the  optically  pumped 
alkali-metal  vapor  (Cs,  Rb  or  K)  in  the  spin-exchange  ceU.  This  work  has  shown  that 
many  diode  laser  arrays  have  a  very  nonuniform  spatial  intensity  distribution.  This 
problem  can  be  solved  by  using  fiber-coupled  arrays,  though  the  laser  cost  is  stiU  an 
issue.  We  are  working  with  the  laser  manufacturer.  Opto  Power,  Inc.,  to  develop  cost- 
efiective  lasers.  More  detail  of  the  new  imaging  method  are  described  in  the  paper 
attached  as  Appendix  B.  This  is  the  work  proposed  in  Item  1.  of  the  Statement  of 
Work,  “Diode  lasers.” 

•  A  cryogenic  accumulator  for  laser-polarized  ^^®Xe  has  been  successfully  developed 
to  produce  the  large  amounts  (liters)  of  ^^®Xe  needed  for  diagnostic  studies  of  the 
human  body.  The  accumtdator  is  very  similar  to  the  device  proposed  in  Item  2.  of 
the  Statement  of  Work,  “Large-volume  source  of  spin  polarized  Xe-129.”  A  patent  on 
the  accumulator,  has  been  issued  to  Princeton  University  and  a  description  has  been 
published  in  the  open  literature® . 

•  We  have  completed  a  quantitative  study  of  image  aberrations  due  to  hindered  diffusion 
of  polarized  gas  near  boundaries.  Some  of  the  results  are  included  in  the  publication 
®  attached  as  Appendix  C. 

•  The  first  ^^®Xe  images  of  a  human  lung  were  obtained  in  collaboration  with  the  Ra¬ 
diology  Department  of  the  University  of  Virginia^®.  This  kind  of  experiment  was 
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proposed  in  Item  3.  of  the  Statement  of  Work,  “Magnetic  resonance  imaging  with 
Xe-129.” 

•  Stimulated  by  our  work,  the  first  international  workshop  on  inert  gas  imaging  was 
held  in  Les  Houches,  France  during  the  fall  of  1996. 
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APPENDIX  A 
Statement  of  Work 

for  ARPA  Grant  DAMD  17-94-J-4469 

Laser  Spin-Exchange  Polarized  ®He  and  ^^®Xe 
for  Diagnostics  of  Gas-Permeable  Media 
with  Nuclear  Magnetic  Resonance  Imaging 

1.  Diode  lasers.  We  will  optimize  the  design  of  spin-exchange  optical  pumping  systems 
for  inexpensive  diode  laser  arrays.  The  important  issues  are: 

(a)  The  large  spectral  linewidth  («  2  nm  FWHM)  of  the  diode  lasers.  This  makes 
it  necessary  to  use  pumping  chambers  with  high  gas  pressures,  typically  sev¬ 
eral  atmospheres  or  more,  to  provide  enough  pressure  broadening  of  the  atomic 
absorption  line  of  the  alkali-metal  atoms  to  utilize  most  of  the  laser  light. 

(b)  We  wiU  assess  the  long-term  reliability  of  diode  laser  arrays  for  spin-exchange  op¬ 
tical  pumping.  For  pumping  the  vapor  of  rubidium  metal,  GaAlAs  diode  material 
is  needed.  In  the  past,  the  tendency  of  the  A1  to  oxidize  at  the  end  facets  and 
the  related  propagation  of  dark  line  defects  have  limited  the  lifetime  of  GaAlAs 
diode  lasers.  However,  major  improvements  in  processing  may  have  resolved  these 
problems.  If  GaAlAs  laser  lifetime  should  be  a  serious  issue,  there  are  good  rea¬ 
sons  to  expect  much  longer  lifetimes  from  aluminum-free  materials  like  GaSbAs, 
which  can  be  grown  by  MO  VCD.  The  operating  wavelengths  of  aluminum-free 
lasers  are  too  long  to  pump  Rb  vapor,  but  they  would  be  an  excellent  match  for 
Cs  vapor. 

2.  Large-volume  source  of  spin  polarized  Xe-129.  So  far,  only  He-3  is  available 
in  sufficient  quantities  (liters)  for  magnetic  resonance  imaging  of  human  lungs.  We 
will  develop  large  volume  sources  of  Xe-129,  which  is  of  interest  because  of  its  much 
greater  solubility  in  human  tissue  than  He-3,  and  because  there  is  an  inexhaustible, 
readily  recoverable  supply  from  the  atmosphere.  The  issues  we  must  address  here  are: 

(a)  The  spin-exchange  and  spin  destruction  cross  sections  of  Xe-129  with  alkah-metal 
atoms  are  some  five  orders  of  magnitude  bigger  than  those  for  He-3.  Other  things 
being  equal,  this  means  that  Xe-129  can  be  polarized  five  orders  of  magnitude 
more  quickly  than  He-3.  Unfortunately,  Xe-129  is  also  five  orders  of  magnitude 
more  potent  in  destroying  the  spin  of  the  alkali-metal  atoms  than  He-3.  Selecting 
the  optimum  concentration  of  Xe-129  and  other  gases  in  the  optical  pumping 
chamber  is  therefore  crucial  to  allow  fast  pumping  while  maintaining  adequate 
spin  polarization  of  the  alkali-metal  atoms. 

(b)  It  is  more  difficult  to  store  spin  polarized  Xe-129  for  long  periods  of  time  than 
is  the  case  for  He-3  gas,  which  can  be  kept  polarized  at  room  temperatures  for 
many  hours.  We  expect  to  be  able  to  solve  this  problem  by  developing  deuterated 
wall  coating  materials  for  the  cells  that  contain  Xe-129.  Work  in  our  laboratory 
has  shown  that  most  of  the  spin  relaxation  of  Xe-129  is  due  to  interactions  with 
the  nuclear  magnetic  moments  of  protons  of  hydrogenated  wall  materials.  Since 
deuterons  have  a  much  smaller  magnetic  moment  than  protons,  they  will  cause 
less  spin  depolarization. 
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(c)  A  very  convenient  way  to  maintain  spin  polarized  Xe-129  for  long  periods  of  time 
is  to  freeze  it  at  liquid  nitrogen  temperatures  or  below,  where  it  has  completely 
reproducible  spin  relaxation  times  of  many  hours.  We  wiU  incorporate  cryogenic 
storage  of  Xe-129  in  our  systems. 

3.  Magnetic  resonance  imaging  with  Xe-129.  Although  Xe-129  was  used  in  the 
very  first  magnetic  resonance  images  with  laser-polarized  gases  -  a  mouse  lung  imaged 
in  a  collaboration  between  our  Princeton  group  and  a  group  at  the  State  University 
of  New  York  at  Stony  Brook  -the  most  spectacular  recent  images  have  been  obtained 
with  He-3,  which  until  now  has  been  available  in  much  larger  volumes  than  Xe-129. 
With  the  success  of  our  efforts  to  make  much  larger  amounts  of  spin-polarized  Xe-129, 
we  will  use  it  to  image  laboratory  animals  and  perhaps  humans.  Xenon  gas  is  very 
soluble  in  blood,  and  there  is  a  good  possibility  that  enough  Xe-129  will  dissolve  in 
the  blood  to  permit  magnetic  resonance  imaging  of  body  tissue.  This  work  will  be 
done  in  collaboration  with  enthusiastic  research  partners  at  the  medical  schools  of 
Duke  University,  the  University  of  Virginia,  the  University  of  Texas  at  San  Antonio, 
and  Vanderbilt  University. 
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Three-dimensional  imaging  of  spin  polarization  of  alkaii-metai  vapor  AFVEmix  b 
in  opticai  pumping  ceiis 

A.  R.  Young,  S.  Appelt,  A.  Ben-Amar  Baranga,  C.  Erickson,  and  W.  Happer 

Department  of  Physics,  Princeton  University,  Princeton,  New  Jersey  08544 
(Received  2  January  1997;  accepted  for  publication  11  April  1997) 

We  present  detailed  images  of  the  Rb  polarization  in  a  high-pressure  optical  pumping  cell.  Images 
taken  at  high  magnetic  field  strengths  ('=^37  G)  yield  precise  information  on  the  relative  occupation 
probability  of  the  hyperfine  Zeeman  sublevels.  Low-field  (^4  G)  images  yield  very  high  resolution 
(less  than  1  mm)  maps  of  the  spatial  distribution  of  the  Rb  polarization.  For  optical  pumping  in  the 
presence  of  high-pressure  buffer  gases,  the  measured  populations  of  the  hyperfine  magnetic 
sublevels  are  shown  to  be  well  described  by  a  spin  temperature,  even  though  the  Rb-Rb  spin 
exchange  rate  is  much  smaller  than  the  rates  of  optical  pumping  or  electron-randomizing  Rb-Xe 
collisions  in  the  buffer  gas.  ©  1997  American  Institute  of  Physics.  [50003-6951(97)02223-7] 


The  high-volume  production  of  laser-polarized  ^^^Xe 
and  ^He  provides  the  basis  for  a  growing  number  of  experi¬ 
mental  and  clinical  techniques.  Medical  applications  of  mag¬ 
netic  resonance  imaging  with  laser-polarized  gases  have 
stimulated  the  development  of  spin-exchange  optical  pump¬ 
ing  systems  capable  of  producing  one  or  two  atm-liters/hour 
of  highly  polarized  ^He  and  ^“^Xe.^  These  high-volume  po¬ 
larization  systems  are  in  a  state  of  rapid  evolution,  and  tools 
are  required  to  precisely  characterize  and  optimize  their  per¬ 
formance.  One  of  the  most  important  properties  of  such  sys¬ 
tems  is  the  three-dimensional  distribution  of  absolute 
electron-spin  polarization  P  of  the  alkali-metal  atoms  in  the 
optical  pumping  cell.  In  this  letter  we  describe  methods 
which  make  such  measurements  possible. 

Although  our  methods  may  be  utilized  over  a  wide  va¬ 
riety  of  experimental  conditions,  here  we  discuss  measure¬ 
ments  of  the  Rb  polarization  in  an  optical  pumping  cell  con¬ 
taining  7.6  atm  (STP)  of  gas,  composed  of  98%  "^He,  1.3% 
Xe  (of  natural  isotopic  composition),  about  0.7%  N2,  and 
saturated  Rb  vapor  in  equilibrium  with  a  few  milligrams  of 
natural  Rb  metal.  The  natural  isotopic  fractions  of  Rb  are 
72%  ^-"Rb  (7=5/2)  and  28%  ®^Rb  (7=3/2),  where  7  is  the 
nuclear  spin.  This  gas  mixture  has  been  used  to  produce 
large  volumes  of  ^^^Xe  in  a  flowing-gas  apparatus.^  The  N2 
prevents  radiation  trapping  by  quenching  the  fluorescence  of 
excited  Rb  atoms. ^  The  "^He  pressure  broadens  the  Lorentz- 
ian  absorption  profile  of  the  Rb  atoms  to  0.27  nm  full  width 
at  half-maximum  (FWHM)  to  better  match  the  2  to  4  nm 
FWHM  emission  profiles  of  the  various  diode  laser  arrays 
used  in  this  work.  The  D1  absorption  profile  for  the  Rb  was 
measured  and  served  as  an  independent  determination  of  the 
"^He  gas  pressure  in  the  cell.  The  high-pressure  "^He  also  acts 
as  a  buffer  gas,  greatly  suppressing  diffusion  of  the  polarized 
Rb  atoms  (the  diffusion  coefficient  at  an  operating  tempera¬ 
ture  of  100  °C  is  about  0.06  cm^/s  for  Rb).  In  such  cells,  the 
Rb  polarization  is  generated  by  optical  pumping  with  circu¬ 
larly  polarized  light  at  a  mean  rate  R.  The  polarization  de¬ 
struction  is  dominated  by  electron-randomizing  collisions 
with  Xe  atoms  at  a  rate  l/Txe.  Other  relaxation  mechanisms 
for  the  Rb  are  relatively  minor. 

For  the  high-pressure  conditions  characteristic  of  spin- 
exchange  optical  pumping,  the  spin  sublevels  \Fm)  of  total 
spin  angular  momentum  F  =  I±II2  and  azimuthal  quantum 
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number  m  are  populated  with  a  probability  very  close  to  the 
spin-temperature  limit: 

piFm)  =  e^^lZ,  (1) 

where  the  Zustandsumme  is  Z  =  The  spin- 

temperature  parameter  is  yS— 2  tanh”^  The  importance 
of  the  spin-temperature  distribution  for  alkali-metal  atoms 
was  first  pointed  out  by  Anderson  et  al?  for  conditions  of 
rapid  spin-exchange  collisions  between  pairs  of  alkali-metal 
atoms.  It  is  a  remarkable  fact,  which  we  will  justify  with 
reference  to  experiments  below,  that  the  spin-temperature 
distribution  also  prevails  in  high-pressure,  spin-exchange  op¬ 
tical  pumping  cells,  irrespective  of  whether  the  rate  of 
Rb-Rb  spin  exchange  collisions. 

The  experimental  apparatus  is  depicted  in  Fig.  1(a).  The 
cell,  a  glass  cylinder,  with  i.d.  of  2.4  cm  and  length  of  8.cm, 
is  placed  in  a  hot-air  oven  to  maintain  adequate  Rb  number 
densities  for  optical  pumping.  The  cell  is  pumped  by  a  diode 
state  laser.  For  most  of  this  work,  we’ used  a  water-cooled 
stack  of  AlGaAs/GaAs  laser-diode  arrays  with  an  overall 
spectral  width  of  roughly  4  nm  FWHM  and  approximately 
80  W  of  laser  power  at  the  795  nm  Rb  D1  line.  The  laser 
light  traverses  two  cylindrical  lenses  (LI  and  L2),  a  polariz¬ 
ing  beam-splitter  cube,  and  a  circular  polarizer.  It  enters  the 
cell  parallel  to  a  uniform  magnetic  field  Bq  which  defines  the 
z  axis  of  a  coordinate  system.  The  x  axis  of  the  coordinate 
system  is  defined  by  a  15  mW,  circularly  polarized,  Ti:Sap- 
phire  laser  probe  beam,  which  enters  the  cell  at  right  angles 
to  Bq  .  The  transmitted  intensity  of  the  probe  laser  is  moni¬ 
tored  by  a  photodiode  operating  in  photoconductive  mode. 
Concentric  with  the  probe  beam  are  a  pair  of  coils  on  either 
side  of  the  cell  in  a  Helmholtz  configuration.  When  radio 
frequency  (RF)  current  is  applied  to  the  coils  and  tuned  to 
the  Bohr  condition  for  transitions  between  the  Zeeman  sub- 
levels  \Fm)  and  \F,m  —  l),  coherence  and  an  associated 
transverse  magnetization,  rotating  in  synchronism  with  the 
RF  [see  Fig.  1(b)],  is  generated.  As  a  label  of  the  resonant 
coherence  of  the  coupled  sublevels  we  use  their  common 
quantum  number  F  and  their  mean  azimuthal  quantum  num¬ 
ber  in  —  m  —  1/2. 

The  absorption  cross  section  for  D1  light  from  the  probe 
laser  is  given  by  (cr)  =  cro(l  “2S' (S)),  where  s  is  the  mean 
photon  spin  of  the  probe  light  and  (S)  is  the  mean  electronic 
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M(a))=MoX 


Fin 

[jfm  sin  (Ot- ((O- (Of ^)cos  <ot] 


(2) 


Here  Mq  is  a  function  of  the  RF  power  and  the  probe  laser 
power,  <0  is  the  RF  frequency,  and  (Of^  is  the  transition 
frequency.  Statistical  weights  of  spin  multiplets  are  denoted 
by  [F]  =  2F+1,  [/]  =  2/+l,  etc.  The  damping  rate  of  the 
coherence  is,  with  .^=(R  +  l/Fxe), 


3[f]^+l-4m^-4/w[/]PRb(-l) 


IfVl-F 


The  population  differences  of  the  sublevels  arc  incorporated 
into  the  factor  Apfm  =  [p(F,»i+ l/2)-p(F,/n- 1/2)], 


2e^"'  sinh^  (l3/2) 

^PFm  gjjjjj  p  j-/]^ 


(4) 


FIG.  1.  (a)  Experimental  apparatus,  as  described  in  the  text.  The  total  mag¬ 
netic  field  B  is  the  resultant  of  a  uniform  field  Bq  and  a  linear  gradient  field 
(r- VB),  both  along  the  z  axis.  Not  depicted  are  the  coils  used  to  generate 
the  gradient  fields  and  the  oven,  (b)  Schematic  sketch  of  imaging  with  a 
linear  gradient.  The  lock-in  signal  arises  solely  from  the  Rb  spins  in  reso¬ 
nance.  As  Bq  is  swept  linearly  in  time,  the  volume  element  in  resonance  is 
swept  along  the  x  axis. 

spin  of  the  Rb.  As  it  rotates  into  and  out  of  alignment  with 
the  probe  laser  axis,  the  rotating  transverse  spin  (S^)  asso¬ 
ciated  with  the  coherence  Fm  modulates  the  transmission  of 
the  probe  laser. 

As  in  nuclear  magnetic  resonance  imaging,'^  information 
about  the  spatial  variation  of  the  polarization  along  the  probe 
beam  can  be  obtained  by  adding  a  spatial  gradient  to  the  field 
Bq  [see  Fig.  1(b)].  This  gradient  is  generated  by  coils  ar¬ 
ranged  to  produce  a  constant  gradient  {dBJdx)  along  the 
probe  beam  axis.  Under  these  conditions,  the  rotating  trans¬ 
verse  magnetization  is  produced  on  thin  surfaces,  perpen¬ 
dicular  to  the  probe  beam.  Sweeping  the  magnitude  of  Bq 
causes  the  surfaces  to  move  along  the  probe  beam  from  one 
side  of  the  cell  to  the  other.  Thus,  the  modulation  amplitude 
is  a  measure  of  the  polarization,  P Rb~  2(^2),  at  the  resonant 
spatial  location  x  along  the  probe  beam.  The  probe  beam 
pierces  theyz  plane  at  the  coordinates  (y,z).  Since  the  probe 
beam  is  only  a  few  mm  or  less  in  diameter,  the  spatial  coor¬ 
dinates  of  the  polarization  P Rb  that  generate  the  modulation 
signal  are  determined  within  a  volume  element  on  the  order 
of  1  mm^. 

In  this  work  we  explored  two  different  imaging  regimes: 
a  “high-field”  regime  G)  where  the  different  Zee- 

man  transitions  are  well  resolved  from  each  other,  and  a 
“low-field”  regime  (Fo'“4  G)  where  all  of  the  Zeeman 
transitions  occur  at  nearly  the  same  frequency.  High-field 
measurements  with  no  field  gradients  provide  detailed  infor¬ 
mation  concerning  the  population  of  the  magnetic  sublevels. 
For  Rb  vapor  in  the  spin  temperature  limit  and  for  suffi¬ 
ciently  weak  RF  field  amplitude,  the  theoretical  modulation 
amplitude  of  the  probe  beam  is^ 


In  this  letter  we  discuss  images  of  ®^Rb  polarization.  Similar 
results  have  been  obtained  for  *''Rb.  The  remarkable  agree¬ 
ment  between  theory  and  experiment,  shown  in  Fig.  2(a), 
demonstrates  that  the  population  distributions  are  indeed 
close  to  the  spin  temperature  limit. 

RF  resonances  with  no  field  gradient  and  uniform  illu¬ 
mination  of  the  vapor  may  be  viewed  as  point-spread  func¬ 
tions  for  gradient  imaging.  The  width  is  dominated  by  the 
optical  pumping  rate  R  and  the  electron  randomization  rate 


37.5  38  B  (Gauss) 


FIG.  2.  (a)  “Rb  signal  in  the  absence  of  an  imaging  gradient.  T^e  signal 
amplitudes  of  the  resonance  peaks  are  in  excellent  agreement  with  a  spin 
temperature  distribution.  Also  depicted  is  a  fit  of  the  data  to  the  functional 
form  of  Eq.  (2).  (b)  “Rb  signal  amplitude  under  similar  conditions,  except 
that  an  imaging  gradient  of  1  G/cm  has  now  been  applied  (only  the  first  two 
resonances  are  presented).  The  ratio  of  these  two  resonance  peaks  can  be 
used  to  make  a  map  of  the  absolute  Rb  polarization,  as  is  depicted  in  the 
inset. 

Young  et  at. 
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FIG.  3.  (a)  The  point  spread  function  obtained  from  a  low-field  magnetic 
sweep  measurement  (no  gradient  present),  indicating  a  spatial  resolution  of 
about  0.7  mm.  (b)  Low-field  images:  in  the  upper  trace  the  cell  was  illumi¬ 
nated  homogeneously,  while  in  the  lower  trace  diverse  peaks  appear  in  the 
image,  reflecting  the  stacked  bar  structure  of  the  diode  laser  imaged  opti¬ 
cally  into  the  cell. 

l/Jxe  ‘  We  have  explored  the  dependence  of  the  resonance 
width  on  R,  VTxq,  and  less  important  variables.  We  find 
excellent  agreement  with  the  theoretical  predictions.^ 

For  imaging  in  high  magnetic  fields  with  resolved  Zee- 
man  resonances,  the  field  gradient  must  not  be  so  big  that 
more  than  one  Zeeman  resonance  is  excited  in  the  cell.  This 
sets  an  upper  limit  on  the  gradient  one  may  add  to  Bq  .  Fig¬ 
ure  2(b)  shows  two  resonances  taken  with  a  gradient  some¬ 
what  less  than  the  limiting  value  and  very  small  illuminated 
volume  (2.5  mm  FWHM  pumping  beam).  The  ratio  of  the 
amplitudes  of  these  resonances,  at  corresponding  points  in 
the  cell,  can  be  converted  into  a  map  of  the  absolute  polar¬ 
ization  of  the  Rb,  with  a  resolution  of  less  than  1.2  mm  along 
the  probe  laser  path  [see  inset.  Fig.  2(b)]. 

Although  reasonable  imaging  can  be  performed  at  high 
field,  much  better  resolution  and  higher  signal-to-noise  ratios 
can  be  achieved  by  reducing  Bq  io  a  value  where  the  differ¬ 
ences  in  the  resonance  frequencies  (the  “quadratic  Zee- 
man  splittings”)  are  small  compared  to  the  damping  rates 
ypffi.  Then  there  will  be  only  one  resonance  frequency  for 
each  Rb  isotope,  and  relatively  large  field  gradients  can  be 
used.^  The  theoretical  analysis  of  the  low-field  images  is 
more  difficult,  because  the  various  Zeeman  resonances  are 
coupled  to  each  other,  but  the  low-field  image  signals  turn 
out  to  be  very  nearly  linear  in  . 

For  example,  in  Fig.  3  we  present  low-field  scans  along 
the  laser  beam  under  a  variety  of  illumination  conditions.  In 
Fig.  3(a)  we  depict  the  polarization  signal  with  no  gradient 
applied.  This  trace  yields  an  estimate  of  0.7  mm  for  our 
spatial  resolution.  In  Fig.  3(b)  we  present  the  polarization 
signal  from  a  uniformly  illuminated  cell,  illustrating  the  rela¬ 
tively  uniform  polarization  distribution  which  results.  Also 
depicted  in  Fig.  3(b)  is  the  polarization  distribution  produced 
when  the  light  from  a  stack  of  laser-diode  arrays  was  simply 


FIG.  4,  Two-dimensional,  low-field  image  for  a  focus  of  the  pumping  beam 
provided  by  a  cylindrical  lens  in  front  of  the  cell. 

imaged  in  the  cell,  with  no  attempt  to  provide  uniform  illu¬ 
mination.  The  extreme  variation  of  the  polarization  is  pro¬ 
duced  by  bands  of  intense  light  coming  from  individual 
laser-diode  arrays.  This  trace  also  demonstrates  that  the  Rb 
polarization  is  “frozen”  in  place  in  high-pressure  optical 
pumping  ceils. 

By  taking  successive  scans  at  low  field,  it  is  possible  to 
produce  two-  and  three-dimensional  images  of  the  polariza¬ 
tion.  In  Fig.  4  we  present  an  image  of  the  polarization  in  the 
xz  plane  under  low  field  conditions.  The  laser  beam  was 
brought  to  a  focus  in  the  cell,  to  emphasize  that  extreme 
spatial  variations  in  the  polarization  can  be  produced  and 
imaged  in  high-pressure  optical  pumping  cells. 

Our  technique  is  applicable  under  a  wide  variety  of  con¬ 
ditions,  and  can  be  generalized  in  a  straightforward  way  to 
produce  images  of  the  absolute  polarization  of  the  Rb  over 
the  entire  volume  of  an  optical  pumping  cell.  The  excellent 
signal-to-noise  ratios  permit  submillimeter  resolution  of  the 
spatial  variations  of  the  Rb  polarization.  This  imaging  tech¬ 
nique  should  provide  a  valuable  tool  in  the  characterization 
of  the  optical  pumping  process  and  in  the  development  of 
more  effective  means  of  polarizing  large  volumes  of  noble 
gas  atoms. 

This  work  was  supported  by  the  AFOSR  and  the  De¬ 
fense  Advanced  Research  Project  Agency,  with  assistance 
from  the  Princeton  cyclotron  group,  funded  by  the  NSF. 
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Abstract 

We  have  produced  one-dimensional  (1-D)  magnetic  resonance  images  of  hyperpolarized  ’He  in  square  glass  cells  and 
demonstrated  edge  enhancement  of  the  signal  intensity  near  the  parallel  impermeable  boundaries.  The  size  and  position 
of  the  edge-enhancement  peaks  approximately  agrees  with  the  theoretical  prediction  based  on  solutions  of  the  1-D  Torrey 
equation  with  boundary  conditions.  Additional  distortion  in  our  images  is  due  to  the  long-range  dipole  fields  produced 
by  the  polarized  spins  and  by  the  bulk  magnetic  susceptibility  of  the  glass.  Our  experiments,  done  in  an  appUed  field  of 
//o  =  31 G,  demonstrate  the  low-field  imaging  capacity  of  hyperpolarized  noble  gases. 


1.  Introduction 

Nuclear  magnetic  resonance  imaging  (MRI)  relies 
on  the  application  of  known  gradients  to  encode  spa¬ 
tial  information  in  the  local  Larmor  precession  fre¬ 
quency  of  the  signal-source  nuclei.  At  resolutions  ap¬ 
proaching  the  characteristic  diffusion  length  of  the  sig¬ 
nal  source,  diffusion  may  affect  this  correspondence 
of  frequency  to  position.  The  most  common  MRI  sig¬ 
nal  source,  particularly  for  biological  imaging,  is  'H 
in  water,  for  which  the  critical  length  scale  is  1  -10  fim 
[  1-3] .  Recently,  laser  techniques  for  hyperpolarizing 
inert-gas  nuclei  [4-7]  have  produced  the  first  MR  im¬ 
ages  using  an  inert  gas  as  the  signal  source  [8-12]. 
For  imaging  of  gas  spaces  such  as  the  lungs,  diffusion 
effects  are  potentially  important  on  much  larger  length 
scales,  because  the  diffusion  coefficient  of  a  typical 
gas  at  atmospheric  pressure  is  about  10*  greater  than 
that  for  water. 

Piitz,  Barsky,  and  Schulten  have  elucidated  the 
effects  on  MRI  of  diffusion  in  restricted  volumes, 
particularly  the  effect  known  as  edge  enhancement. 


whereby  the  image  intensity  is  increased  near  the 
impermeable  boundaries  of  an  imaged  region  [13]. 
They  and  co-workers  have  also  observed  the  effect 
in  high-resolution  MR  images  of  water  in  a  thin 
glass  capillary  [14,15].  Hyslop  and  Lauterbur  have 
also  treated  the  problem  theoretically  using  computer 
simulations  [16].  We  present  here  MR  images  of 
hyperpolarized  ’He  in  a  geometry  approximating  a 
one-dimensional  (1-D)  region  with  fixed,  parallel, 
impermeable  boundaries.  MRI  with  hyperpolarized 
gases  offers  a  more  convenient  way  to  study  edge 
enhancement  for  several  reasons.  First,  the  hyperpo¬ 
larized  gases  are  an  intense  signal  source,  typically  at 
least  ten  times  stronger  than  water  in  MRI  fields  [9] . 
Second,  the  more  rapid  gas  diffusion  means  that  the 
range  of  edge  enhancement  extends  farther  from  the 
boundary,  and  the  effect  may  thus  be  observed  with 
smaller  applied  gradients  and  better  signal-to-noise 
ratio  (SNR).  Finally,  the  diffusion  constant  is  easily 
varied  over  a  few  orders  of  magnitude  by  varying  the 
gas  pressure,  making  it  possible  to  observe  the  effect 
over  a  broad  range  of  parameters. 
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2.  Theory 


We  consider  a  1-D  imaging  experiment  along  the 
2 -axis,  which  begins  by  uniformly  tipping  the  nuclear 
spins  in  a  uniformly  magnetized  sample  region  by  an 
angle  6  with  a  resonant  rf  pulse.  At  time  ?  =  0,  the  uni¬ 
form  rotating-frame  transverse  magnetization,  = 
Mx  *f  iA/y,  is  given  by 

M+(z,0)  -  Pnylhsind,  (1) 


where  P  is  the  polarization,  n  is  the  density  of  nu¬ 
clei,  /  is  the  nuclear  spin,  and  y  is  the  gyromagnetic 
ratio.  For  t  >  0,  the  spins  precess  freely  in  the  plane 
transverse  to  the  uniform  applied  field  Hq  =  Hqz,  in 
the  presence  of  a  single  constantly  applied  gradient 
Gz  =  A  1-D  MR  image  is  completely  repre¬ 

sented  by  the  Fourier  spectrum  of  the  subsequent  free 
induction  decay  (FID).  Conceptually,  this  is  perhaps 
the  simplest  MRI  experiment,  and  it  is  the  best  model 
for  the  work  presented  here.  Moreover,  it  has  a  partic¬ 
ular  relevance  to  gas-phase  MRI.  In  MRI,  it  is  ad¬ 
vantageous  to  acquire  an  image  from  one  or  more  spin 
echoes.  In  gas-phase  MRI,  however,  echo  techniques 
are  more  difficult  to  implement,  because  rapid  diffu¬ 
sion  severely  attenuates  the  echo  intensity.  Techniques 
such  as  projection  acquisition  [17,18],  which  can  ac¬ 
quire  an  image  directly  from  the  FID,  may  eventually 
prove  more  useful  for  gas-phase  MRI. 

The  evolution  of  Af  +  for  r  >  0  was  first  treated  by 
Torrey  [19].  In  the  1-D  case  for  a  reference  frame 
rotating  near  the  Larmor  frequency  vo  =  yHo/lTr, 


dM^ 

dt 


Mj, 

=  iyGzZAf+ 

h 


d'^M^ 
dz^  ’ 


(2) 


where  D  is  the  diffusion  coefficient,  and  T2  is  the 
transverse  relaxation  time.  For  an  ideal  gas,  =  7i . 
For  hyperpolarized  inert  gases,  T2  is  usually  so  long 
that  the  corresponding  term  in  (2)  is  negligible.  The 
well-known  solution  to  (2)  in  the  case  of  no  bound¬ 
aries  (unrestricted  diffusion)  was  first  derived  by 
Hahn  [20]  and  is  the  basis  for  measuring  D  in  vari¬ 
ous  liquids.  Stoller,  Happer,  and  Dyson  [21]  included 
boundary  conditions,  treating  the  more  general  case 
of  restricted  diffusion  in  one  dimension.  Drawing 
on  their  results,  deSwiet  and  Sen  [22]  generalized 
to  higher  dimensions  and  provided  a  theoretical  in¬ 
terpretation  of  edge  enhancement,  exploring  various 


regimes  in  which  the  Hahn-echo-amplitude  formal¬ 
ism  breaks  down.  Piitz  and  co-workers  developed  a 
computer  algorithm  to  generate  images  based  on  the 
formal  solution  of  (2)  in  one,  two,  and  three  dimen¬ 
sions  (the  line  shape  function).  Fig.  3  of  Ref.  [13] 
displays  their  computer-generated  1-D  MR  images 
for  various  values  of  D,  to  which  our  data  presented 
below  may  be  compared. 

As  pointed  out  by  deSwiet  and  Sen  [22],  there  are 
two  regimes  of  interest  relative  to  the  characteristic 
sample  dimension  L  and  the  natural  length  scale  a 
derived  from  (2): 


For  L  ^  a,  only  spins  near  local  boundaries  are  af¬ 
fected.  Restricted  diffusion  causes  a  local  increase  in 
the  dephasing  time  T2,  leading  to  a  peak  in  the  image 
near  the  boundary.  In  our  analytical  solution  to  the  1-D 
case  with  one  boundary  (to  be  published) ,  the  peak  is 
located  a  distance  1.04Qf  from  the  wall,  whereas  away 
from  the  boundaries  the  image  appears  as  it  would  for 
D  =  0.  For  L  a,  most  of  the  spins  traverse  the 
sample  many  times  during  image  acquisition.  The  re¬ 
sulting  image  corresponds  to  a  broad  averaging  of 
over  the  entire  region.  In  the  1-D  case,  the  spectrum 
appears  as  a  single  central  peak  which  narrows  fur¬ 
ther  for  increasing  a.  For  image  acquisition  from  the 
FID,  both  in  and  between  these  two  regimes,  the  to¬ 
tal  signal  generated  by  the  spins  (the  integrated  image 
intensity)  is  conserved,  given  the  same  Af+(z,0)  for 
each  image.  Signal  is  not  lost  but  is  redistributed  to 
other  parts  of  the  frequency  spectrum  according  to  the 
particular  solution  of  (2).  This  is  not  generally  true 
for  imaging  with  echoes. 


3.  Experiment 

Box-shaped  cells  (Fig.  1)  made  of  Coming  7740 
(Pyrex™)  glass  were  prepared  on  a  high-purity  gas 
handling  system.  The  cells  are  made  from  square  tub¬ 
ing  with  internal  dimensions  6mmxl0mm,  closed 
with  a  flat  bottom  at  one  end  and  with  a  roughly  coni¬ 
cal  pull-off  at  the  other.  This  additional  volume  at  the 
top  center  of  the  cell  is  an  unavoidable  consequence  of 
sealing  it  away  from  a  glass  manifold  during  its  pro¬ 
duction.  The  cells  contain  10- 100  mg  of  distilled  Rb 
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Fig.  1.  The  two  cells  imaged  in  this  work.  The  main  field  and 
gradient  are  applied  in  the  plane  of  the  photo  from  left  to  right. 
The  depth  2a  =  6  mm,  the  height  2/»  =  12  mm  (not  including  the 
pull-off),  and  the  width  2L  =  1 0mm.  The  cell  on  the  left  con¬ 
tains  2.5  atm  -’He  and  0.3  atm  N2,  yielding  a  diffusion  constant 
D  =  0.59cm^/s  [23,24].  The  cell  on  the  right  contains  463TorT 
He  and  30Toit  N2,  yielding  D  =  2.6cm^/s.  Pressures  are  mea¬ 
sured  at  295  K. 

metal,  1 0-200  Torr  N2  gas,  and  0.2-2.5atm  of  ^He, 
where  the  pressures  are  measured  at  295  K.  Precise 
values  of  the  contents  of  the  two  cells  used  to  provide 
data  for  this  work  are  given  in  Fig.  1.  We  measured  a 
^He  longitudinal  relaxation  time  T\  >  80  hr  for  these 
cells. 

Images  are  generated  by  a  single-coil  single-channel 
pulse-NMR  spectrometer  operating  at  vq  =  101  kHz, 
the  Larmor  frequency  of  ^He  for  /fo  =  3 1 G.  The  main 
field  is  provided  by  Helmholtz  coils  ( 84  cm  dia. ) ,  and 
the  linear  gradient  is  provided  by  a  second  con¬ 
centric  set  of  counter-wound  coils  ( 14.5  cm  dia.)  with 
a  separation-to-radius  ratio  of  \/l  to  suppress  higher 
order  gradients  [25],  The  probe  coil  is  a  solenoid 
(2.5  cm  dia.  by  5.0cm  long)  of  274  turns  wound  in 
two  layers  and  has  an  unloaded  quality  factor  Qo  sa 
40.  The  cells  are  mounted  at  the  center  of  the  probe 
such  that  the  1-D  imaged  region  along  the  z-axis  cor¬ 
responds  to  the  distance  2L  =  10  mm  between  one 
set  of  parallel  cell  walls.  The  spectrometer  transmits  a 
single  manually-triggered  square  rf  excitation  pulse  to 
the  probe  and  receives  the  subsequent  FID.  The  FID 
is  digitized  with  an  oscilloscope  and  downloaded  to  a 
computer,  where  a  Fast  Fourier  Transform  (FFT)  is 
performed  to  obtain  the  frequency  spectrum. 

In  order  that  the  spectrum  closely  represent  a  1-D 
MR  image  of  the  cell,  several  conditions  are  met.  The 


gradient  Gj  is  on  continuously  to  avoid  any  prob¬ 
lems  with  transients,  and  the  image  is  acquired  directly 
from  the  FID.  Both  rf-pulse  duration  and  total  receiver 
deadtime  are  thus  much  smaller  than  where 

Ap  =  (y/7r)GjL  is  the  spread  in  resonance  frequen¬ 
cies  across  the  cell  with  the  gradient  on.  At  these  low 
frequencies,  ringing  of  the  probe  is  especially  a  prob¬ 
lem.  A  resistor  in  parallel  with  the  tuned  tank  circuit 
suppresses  the  ringdown  at  the  expense  of  sensitivity, 
reducing  Q  from  40  to  about  10.  The  combination  of 
low  filling  factor  and  low  Q  reduces  the  sensitivity  by 
more  than  an  order  of  magnitude  over  what  it  would  be 
for  a  well-filled  unloaded  tank  circuit.  Nevertheless, 
even  for  small  tip  angles,  we  still  have  SNR  s»  100 
for  the  2.5  atm  cell  over  a  bandwidth  ^.v  k  100  Hz 
with  the  cell  polarized  to  about  15%.  Finally,  the  rf- 
excitation  field  H\  must  be  uniform  throughout  the 
cell;  we  thus  constructed  an  oversized  coil  for  which 
the  filling  factor  for  the  cells  is  about  2%.  Numerical 
calculations  indicate  that  the  spread  in  H\  across  the 
cell  is  a  few  percent  under  these  conditions.  We  note 
that  the  2.5  atm  cell  must  be  handled  carefully  when 
its  ^He  spins  are  in  the  high-energy  Zeeman  state:  even 
with  the  g-spoiling  resistor  and  low  filling  factor,  if 
the  sample  is  placed  in  the  coil  with  =  0,  a  very 
rapid  and  spontaneous  loss  of  polarization  (masing) 
occurs  [26]. 

The  cell  is  first  polarized  through  optical  pumping 
and  spin  exchange  to  about  15%  in  the  case  of  the 
2.5  atm  cell.  (The  polarization  is  determined  via  a 
comparison  to  the  thermal-equilibrium  proton  signal 
from  water  contained  in  a  similarly  shaped  cell.)  The 
polarized  cell  is  then  removed  to  the  spectrometer.  A 
series  of  images  is  taken  using  low-angle  pulses  at 
various  gradient  strengths,  both  with  the  longitudinal 
magnetization  parallel  ( || )  and  antiparallel  ( | )  to 
Hq.  The  reversal  of  M^z  is  accomplished  periodically 
with  a  nearly  lossless  adiabatic  inversion.  The  series 
is  repeated  several  times  for  successively  lower  polar¬ 
ization.  When  the  polarization  is  nearly  exhausted,  a 
few  large-angle  pulses  are  used  to  improve  the  SNR. 

4.  Results  and  discussion 

Fig.  2a  shows  a  series  of  images  of  the  2.5  atm 
cell  with  an  initial  polarization  of  16%.  In  all  images, 

Gj  >  0  and  higher  precession  frequencies  correspond 
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G^=4.0  mG/mm 


position  (mm) 


Rg.  2.  (a)  Images  of  the  2.5  atm  cell  showing  edge  enhancement  at  three  gradient  strengths  for  the  spin  magnetization  parallel  (left 
column)  and  antiparallel  (right  column)  to  the  applied  field.  The  spectra  in  one  column  were  inverted  relative  to  the  other  for  easier 
comparison  of  features.  The  central  peak  is  due  to  the  pull-off.  The  applied  field  and  applied  gradient  point  to  the  right  in  all  images.  The 
dashed  lines  mark  the  cell  walls.  The  intensity  units  are  arbitrary  but  identical  for  all  values  of  G..  The  polarization  was  about  16%  and 
was  approximately  conserved  for  all  images  by  using  a  low-angle  excitation  pulse,  (b)  Distance  of  edge-enhancement  peaks  in  (a)  from 
the  cell  wall  as  a  function  of  applied  gradient  for  the  2.5  atm  cell.  The  line  corresponds  to  the  theoretical  dependence  1 .04a  in  the  limit 
L'^  a  with  D  =  0.59cm2/s.  Solid  (open)  points  correspond  to  the  magnetization  parallel  (antiparallel)  to  the  applied  field. 


to  the  positive  z-axis.  The  left  (right)  column  of  im^ 
ages  corresponds  to  M^z  ||  (||)  Hq,  For  the  low  flip 
angles  used,  M+(z,0)  is  approximately  constant  for 
all  the  images  in  a  column.  (There  is  a  small  differ¬ 
ence  between  the  columns  due  to  losses  from  the  in¬ 
tervening  adiabatic  reversal.)  To  within  a  few  percent 
experimental  error,  we  find  that  the  integrated  intensity 
of  the  spectra  in  each  column  is  independent  of 
In  addition  to  the  edge-enhancement  peaks,  there  is  a 
central  peak  due  to  the  pull-off  volume.  Fig.  2b  shows 
the  displacement  of  the  edge-enhancement  peaks  from 
the  walls,  which  is  generally  greater  than  the  calcu¬ 
lated  value  of  \.04a  from  (3)  by  about  10%.  The  dis¬ 
crepancy  is  worst  for  small  G^ ,  where  the  condition 
L  »  a  is  not  well  satisfied,  and  may  also  have  to  do 
with  the  image  distortions  described  below. 

The  images  are  generally  distorted,  most  obviously 
in  that  the  peak  heights  are  asymmetric.  The  asymme¬ 
try  reverses  sign  when  M^z  is  reversed  with  respect  to 
Ho,  and  diminishes  with  increasing  G^.  As  shown  in 


Fig.  3,  the  asymmetry  is  also  proportional  to  polariza¬ 
tion  at  fixed  G^.  Even  at  low  polarizations,  there  is  a 
residual  asymmetry  (the  y-intercept  of  the  plot  in  Fig. 
3)  having  the  same  sign  as  for  the  case  of  M^z  ||  ^o- 
Through  various  combinations  of  reversing  the  direc¬ 
tion  of  Gz  and  rotating  the  cell  by  180°  with  respect  to 
the  z-axis,  we  verified  that  the  spectral  asymmetry  has 
little  to  do  with  asymmetries  either  in  the  cell  shape 
or  in  the  applied  gradient  field.  Rather,  the  experimen¬ 
tal  facts  are  consistent  with  the  distortion  arising  from 
the  internal  long-range  field  produced  at  the  location 
of  each  spin  by  all  of  the  other  spins  in  the  sample. 
Known  as  the  dipolar  demagnetizing  field  [27-29] ,  it 
is  both  nonzero  and  nonuniform  across  all  finite  uni¬ 
formly  magnetized  volumes  other  than  ellipsoids. 

We  consider  R^(z),  the  z-component  of  the  dipo¬ 
lar  demagnetizing  field  averaged  over  the  cross  sec¬ 
tional  slice  of  the  cell  at  each  value  of  z.  For  volumes 
possessing  the  symmetry  of  our  samples  (ignoring  the 
pull-off),  B^(z)  is  evenly  symmetric  about  z  =  0. 


18 


B.  Saam  et  al./ Chemical  Physics  Letters  263  (1996)  481-487 


-10  -50  5  10  -10  -50  5  10 

Position  (mm)  Position  (mm) 


Fig.  3.  Images  of  the  2.5  atm  cell  with  the  same  applied  gradient 
Gz  =  4.0  mG/ mm  but  different  polarizations.  The  magnetization  is 
parallel  to  the  applied  field  for  all  three  images.  The  intensity  units 
are  different  for  each  image  because  larger  excitation  pulses  were 
used  to  increase  SNR  for  the  lower  polarizations.  The  asymmetry, 
measured  as  the  difference  in  edge-enhancement  peak  heights 
divided  by  their  sum,  is  plotted  versus  polarization. 

The  gradient  component  (z)  possesses  a  cor¬ 
responding  odd  symmetry.  In  the  limit  of  zero  diffu¬ 
sion,  the  signal  intensity  /(z),  otherwise  dependent 
only  on  Gj,  is  altered  by  an  amount 

A/(z) 

assuming  ^  and  ignoring  any  dc  shift  of 
the  entire  image  due  to  B^.  Since,  for  even-order  gra¬ 
dients,  Vz5j(-z)  =  -Vzfij(z),  the  signal  lost  at 
some  point  z  =  ^,  for  which  V^B,  (O  and  G,  have 
the  same  sign,  effectively  shows  up  at  z  =  where 
the  two  gradients  have  opposite  signs. 

The  residual  asymmetry,  which  is  observed  directly 
at  low  polarizations,  is  due  to  the  bulk  magnetic  sus¬ 
ceptibility  (BMS)  of  the  glass.  This  asymmetry  is  also 
apparent  at  higher  polarizations  in  that  the  opposing 
images  in  the  two  columns  of  Fig.  2a  are  not  symmet¬ 
ric  with  respect  to  each  other.  The  glass  BMS  gives 
rise  to  a  field  ,  analogous  to  but  independent 
of  the  magnitude  and  orientation  of  We  inde¬ 
pendently  verified  our  sensitivity  to  Bf  by  comparing 
spectra  for  =  0  of  the  2.5  atm  cell  with  a  similar 
(slightly  larger)  spherical  cell  having  about  the  same 


-4  -2  0  2  4 

cell  position  (mm) 


Rg.  4.  The  function  Sf  (z),  the  slice-averaged  dipolar  demagne¬ 
tizing  held,  for  spin  magnetization  parallel  (negative  curvature) 
and  antiparallel  (positive  curvature)  to  the  applied  field,  obtained 
by  numerical  integration  of  (5).  The  applied  field  and  applied 
gradient  point  to  the  right.  Also  shown  is  the  function  fi?(z)  due 
to  the  glass  BMS. 

pressure  of  ^He.  For  an  ideal  sphere,  there  are  no  in¬ 
ternal  fields  due  to  glass  BMS.  The  spectra  for  both 
cells  were  taken  with  large  flip  angles  and  P  <  1%, 
so  the  effects  of  Bf  were  negligible.  The  field  homo¬ 
geneity,  as  given  by  the  transverse  dephasing  time 
was  more  than  a  factor  of  two  better  for  the  spherical 
cell  than  for  the  box-shaped  cell. 

Fig.  4  shows  the  results  of  numerical  calculations 
of  Bf  (z)  for  the  2.5  atm  cell  with  M^z  ||  Hg  and 
16%  polarization.  Also  shown  is  a  similar  curve  for 
the  internal  field  Bf(z)  due  to  the  glass  BMS.  The 
calculations  are  based  on  solutions  of  Laplace’s  equa¬ 
tion  for  the  internal  ff-field  generated  by  the  uniformly 
magnetized  gas  and  glass,  with  respective  surface- 
charge  densities  appropriately  ascribed  [30].  Relative 
to  ffg  >  0,  the  field  B^  is  found  by  numerically  inte¬ 
grating 
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Fig.  5.  An  image  of  the  463  Ton*  cell  with  applied  gradient 
Gz  =  1.5mG/mm.  The  polarization  was  less  than  5%,  and  a 
laige-angle  {0  >  45®)  excitation  pulse  was  used. 

where  2a  and  2b  are  the  transverse  dimensions  of  the 
cell.  A  similar  expression,  in  which  there  is  no  47r/3 
term  and  we  account  for  the  glass  thickness,  yields 
where  and  the  susceptibility  x  = 

- 1 . 14x  10“*^  for  SiOa  [31  ] .  In  Fig.  4,  the  larger  slopes 
near  the  edges  represent  the  gradients  which  give  rise 
to  the  observed  peak  asymmetries.  The  curvatures  of 
B^(z)  and  Bf(z)  correctly  predict  the  sign  of  the 
asymmetry  for  all  images.  In  particular,  the  two  curves 
add  constructively  for  M^z  ||  and  partially  cancel 
for  MzZ  II  ^o»  corresponding  to  a  larger  asymmetry  in 
the  left  column  of  Fig.  2a  than  for  the  corresponding 
( same  Gz )  spectra  in  the  right  column.  The  magnitude 
of  the  asymmetry  depends  on  the  size  of  Vz  ) 

relative  to  Gz .  Hence  the  asymmetry  is  diminished  for 
larger  Gz  at  fixed  polarization  (Fig.  2)  or  for  smaller 
polarization  at  fixed  Gz  (Fig.  3). 

Using  (4)  and  the  slopes  in  Fig.  4,  one  may  also 
estimate  the  magnitude  of  the  peak  asymmetry.  This 
estimate  relies  on  the  condition 

+  (6) 

in  which  case  the  qualitative  features  of  the  image  are 
assumed  unaltered,  while  the  intensity  at  each  point 
changes  according  to  (4) .  We  find  that  the  peak  asym¬ 
metries  estimated  in  this  fashion  are  typically  a  factor 
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of  three  smaller  than  those  measured  directly  from  our 
data.  At  least  some  of  this  discrepancy  is  due  to  the 
fact  that  we  have  ignored  the  complicated  effects  of 
B^  “f  Bf  on  the  particular  solution  to  (2).  Regard¬ 
less,  the  dependence  of  the  sign  and  magnitude  of  the 
asymmetry  on  the  sign  and  magnitude  of  MzZ  leave 
little  doubt  that  the  field  Rf  +  Bf  is  the  principal 
source  of  image  distortion  after  edge  enhancement. 

The  greatest  uncertainty  in  characterizing  our  im¬ 
ages  results  from  our  inability  to  locate  independently 
the  geometric  center  of  the  sample  in  frequency 
space.  Ambient  field  inhomogeneities  and  fluctua¬ 
tions  (about  a  part  in  10“"^ ),  small  dc  components  in 
both  the  applied  and  internal  gradients,  and  the  com¬ 
plicated  nature  of  the  image  distortion  near  the  walls 
all  combine  to  limit  such  a  determination  from  the 
spectra  themselves.  We  chose  to  define  the  center  as 
the  midpoint  between  the  edge-enhancement  peaks, 
recognizing  that  B^  -\-Bf,  in  addition  to  changing  the 
peak  heights,  may  cause  the  peaks  to  shift  asymmet¬ 
rically.  This  is  especially  true  for  small  Gz ,  where  the 
condition  (6)  is  perhaps  only  marginally  satisfied. 

Fig.  5  shows  a  1-D  image  of  the  463-Torr  cell,  which 
has  a  diffusion  coefficient  of  2.6cm^/s  for  ^He.  The 
edge-enhancement  peaks  in  these  images  have  moved 
more  than  halfway  toward  the  center  of  the  cell  but 
have  not  yet  merged  into  a  single  peak.  The  images 
represent  the  transition  region  a  between  the  ex¬ 
treme  conditions  discussed  in  Section  2.  Cells  with 
even  lower  pressures  were  prepared  to  explore  the  the 
region  L  <  a,  but  we  had  technical  problems  produc¬ 
ing  sufficient  polarization  in  them  to  get  good  images. 

We  have  demonstrated  that  1-D  MR  images  of  hy- 
perpolarized  noble  gases  readily  exhibit  edge  enhance¬ 
ment.  Good  images  are  obtained,  even  at  low  Hq,  be¬ 
cause  the  laser-produced  hyperpolarization  is  indepen¬ 
dent  of  Hq.  In  addition,  we  have  observed  and  charac¬ 
terized  two  additional  sources  of  image  distortion:  the 
dipolar  demagnetizing  field  B^  due  to  the  hyperpo- 
larized  spins  and  the  internal  field  Bf  due  to  the  BMS 
of  the  glass  container.  All  of  these  effects  are  present 
or  have  analogs  in  biological  imaging  of  noble  gases. 
Edge  enhancement  or  some  form  of  gradient  averag¬ 
ing  could  show  up  in  the  smaller  airways  of  the  lung. 
The  effects  of  lung-tissue  BMS  on  alveolar  *^^Xe  gas 
in  mice  have  already  been  observed  [32].  At  higher 
^He  polarizations,  the  effect  of  the  dipolar  demagne¬ 
tizing  field  in  a  nonspherical  compartment  might  be 
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observed  even  at  pressures  below  1  atm. 


Acknowledgements 

We  thank  W.S.  Warren  and  d.A.  Johnson  for  help¬ 
ful  discussions,  and  CJ.  Erickson  for  help  in  design¬ 
ing  and  constructing  the  experimental  apparatus.  This 
work  was  supported  by  AFOSR  grant  No.  F49620-94- 
1-0466  and  AREA  grant  No.  DAMD17-94-J-4469. 

References 

[  1  ]  RT.  Callaghan  and  C.D.  Eccles.  J.  Magn.  Res.  78  ( 1988)  I . 

[2]  CB.  Ahn  and  Z.H.  Cho,  Med.  Phys.  16  (1988)  22. 

[3]  C.B.  Ahn  and  Z.H.  Cho,  Magn.  Reson.  Med.  19  (1991)  228. 

[4]  N.D.  Bhaskar  and  W.  Mapper,  Phys.  Rev.  Lett.  49  ( 1982)  25. 

[5]  W.  Mapper,  E.  Miron,  S.  Schaefer,  D.  Schreiber,  W.A.  van 
Wijngaarden  and  X.  Zeng,  Phys.  Rev.  A  29  (1984)  3092. 

[6]  G.D.  Cates,  R.J.  Fitzgerald,  A.S.  Barton,  P.  Bogorad,  M. 
Gatzke,  N.R.  Newbury  and  B.  Saam,  Phys.  Rev.  A  45  (1992) 
4631. 

[7]  T.E.  Chupp,  M.E.  Wagshul,  K.R  Coulter,  A.B.  McDonald 
and  W.  Mapper,  Phys.  Rev.  C  36  (1987)  2244. 

[8]  M.S.  Albert,  G.D.  Cates,  B.  Driehuys,  W.  Mapper,  B.  Saam, 
C.S.  Springer  Jr.  and  A.  Wishnia,  Nature  370  (1994)  199. 

[9]  M.  Middleton,  R.D.  Black,  B.  Saam,  G.D.  Cates,  G.R  Cofer, 
B.  Guenther,  W.  Mapper,  L.W.  Hedlund,  G.A.  Johnson,  K. 
Juwan  and  J.  Swartz,  Magn.  Reson.  Med.  33  (1995)  271. 

[10]  R.D.  Black,  M.  Middleton,  L.W.  Hedlund,  M.D.  Shattuck, 
G.A.  Johnson,  J.  Swartz,  B.  Driehuys,  G.D.  Cates,  G.R  Cofer 
and  W.  Mapper,  Radiology,  submitted  for  publication  ( 1996) . 

[11]  J.R.  MacFall,  M.C.  Charles,  R.D.  Black,  M.  Middleton,  J. 
Swartz,  B.  Saam,  B.  Driehuys,  C.  Erickson,  W.  Mapper,  G.D. 
Cates,  G.A.  Johnson  and  C.E.  Ravin,  Radiology,  submitted 
for  publication  ( 1 996 ) . 


[12]  M.  Ebert,  T.  GroGmann,  W.  Meil,  E.  Otten,  R.  Surkau,  M. 
Leduc,  P.  Bachert,  N.V.  Knopp,  R.  Schad  and  M.  Thelen, 
Lancet  347  (1996)  1297. 

[  13]  B.  Putz,  D.  Barsky  and  K.  Schulten,  J.  Mag.  Res.  97  ( 1992) 
27. 

[14]  B.  Piitz,  D.  Barsky  and  K.  Schulten,  Chem.  Phys.  Lett.  183 
(1991)  391. 

[15]  D.  Barsky,  B.  Piitz,  K.  Schulten,  J.  Schoeniger,  E.W.  Msu 
and  S.  Blackband,  Chem.  Phys.  Lett.  200  (1992)  88. 

[  16]  W.B.  Myslop  and  P.  Lauterbur,  J.  Mag.  Res,  94  (1991)  501. 

[17]  C.B.  Ahn,  S.Y.  Lee,  O.  Nalcioglu  and  Z.M.  Cho,  Med.  Phys. 
13  (1986)  789. 

[18]  S.L.  Gewalt,  G.M.  Glover,  J.R.  MacFall,  L.W.  Medlund  and 
G.A.  Johnson,  Magn.  Reson.  Med.  29  (1993)  99. 

[19]  M.C.  Torrey,  Phys.  Rev.  104  (1956)  563. 

[20]  E.L.  Mahn,  Phys.  Rev.  80  (1950)  580. 

[21]  S.D.  Stoller,  W.  Mapper  and  F.J.  Dyson,  Phys.  Rev.  A  44 
(1991)  7459. 

[22]  T.M.  deSwiet  and  RN.  Sen,  J.  Chem.  Phys.  100  (1994) 
5597. 

[23]  R.  Barbe,  M.  Leduc  and  F.  Laloe,  J.  Phys.  (Paris)  35  ( 1974) 
935. 

[24]  J.C.  Liner  and  S.  Weissman,  J.  Chem.  Phys.  56  ( 1972)  2288. 

[25]  W.A.  Anderson,  Rev,  Sci.  Instr.  32  (1961)  241. 

[26]  T.E.  Chupp,  R.J.  Moare,  R.L.  Walsworth  and  B.  Wu,  Phys. 
Rev.  Lett.  72  (1994)  23. 

[27]  G.  Deville,  M.  Bernier  and  J.M.  Delrieux,  Phys.  Rev.  B  19 
(1979)  5666. 

[28]  R.  Bowtell,  R.M.  Bowley  and  P.  Glover,  J.  Magn.  Reson. 
88  (1990)  643. 

[29]  W.S.  Warren,  S.  Lee,  W.  Richter  and  S.  Vathyam,  Chem. 
Phys.  Lett.  247  (1995)  207. 

[30]  J.D.  Jackson,  Classical  electrodynamics  (Wiley,  New  York, 
1975)  p.  192. 

[31]  D.R.  Lide,  ed.,  Mandbook  of  chemistry  and  physics  (CRC 
Press,  Boca  Raton,  FL,  1992). 

[32]  M.E.  Wagshul,  T.M.  Button,  M.F.  Li,  Z.  Liang,  C.S.  Springer, 
K.  Zhong  and  A.  Wishnia,  Magn.  Reson.  Med.  36  (1996) 
183. 


